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The physics of high harmonics has led to the generation of 
attosecond pulses [lj and to trains of attosecond pulses |2j. Mea- 
surements that confirm the pulse duration are all performed in 
the far field. All pulse duration measurements tacitly assume that 
both the beam's wavefront and intensity profile are independent 
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of frequency. However, if one or both are frequency dependent, 
then the retrieved pulse duration depends on the location where 
the measurement is made. We measure that each harmonic is 
very close to a Gaussian, but we also find that both the intensity 
profile and the beam wavefront depend significantly on the har- 
monic order. Thus, our findings mean that the pulse duration will 
depend on where the pulse is observed. Measurement of spec- 
trally resolved wavefronts along with temporal characterization 
at one single point in the beam would enable complete space-time 
reconstruction of attosecond pulses. Future attosecond science ex- 
periments need not be restricted to spatially averaged observables. 

High harmonic generation is a highly coherent process in which many 
photons are converted into one, leaving the generating atoms in the nonlin- 
ear medium unchanged. High harmonic generation is often described with 
a semiclassical model in which an electron tunnel ionizes, follows a semi- 
classical trajectory and then recombines by emitting a single photon. In 
a gas of many atoms a phase matched output results in laser-like coher- 
ence of the harmonics, which has been confirmed by interfering spectrally 
resolved replica beams [3 J or different spatial regions of the same beam [3]. 
These measurements however, do not address differences in the amplitude 
or curvature of wavefronts between harmonics, even though we know that 
the wavefront depends upon the trajectory [F\. Such a differences will 
be important for all attosecond pulse and pulse train applications. They are 
the subject of this paper. 

To measure the spectrally and spatially resolved amplitudes and phases, 
we used frequency resolved wavefront characterization (or SWORD-Spectral 
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Wavefront Optical Reconstruction by Diffraction) [8]. The measurement is 
free of approximations - it relies only on linear optics. 

Our measurements reveal that the short-trajectory harmonics have an 
essentially Gaussian spatial profile and that the root-mean-square deviations 
from a parabolic wavefront are only < 0.2 radians. They also reveal that 
different harmonics have a significantly different wavefront curvature upon 
emerging from the generating region. Since we measured both amplitude 
and phase of each harmonic, we are able to propagate the fields and deter- 
mine their spatial properties anywhere along the direction of propagation. 
Specifically we determine their waist position (position of the flat phase) 
and size - the virtual origin of the coherent beams. We find that this virtual 
origin of each harmonic is at significantly different position. Particularly im- 
portant, the distance between the waist positions of some harmonics can be 
larger than their respective Rayleigh ranges. Since the harmonics are gen- 
erated in the gas jet, this observation means that they are produced with 
different radii-of-curvature, depending on the harmonic order. It also means 
that, if the harmonics are focused to measure their temporal structure, the 
duration of the pulse will depend on where the measurement is made. 

In our experiment we measure high harmonics produced by focusing an 
800nm, 35fsec, 400/zJ, laser beam onto a Nitrogen gas jet. For SWORD |8j 
a small 20fj, wide slit is scanned across the XUV beam and the spectrally 
resolved diffraction pattern is recorded. Using this diffractogram, we recon- 
struct spectral wavefront phase and amplitudes of each harmonic. As the 
fundamental beam size in the interaction region is much smaller than the 
gas jet size, we can assume a uniform gas density across the beam in the 
generation region and therefore rotational symmetry of the generated har- 
monic beam. The further details about the experimental setup and phase 
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reconstruction are described in the Methods section. Figure [T] (a) shows 
examples of the reconstruction for harmonics 15 and 21. 

The intensity distribution that we measure for different harmonics in Fig. 
[T](a) is Gaussian, within the ~ 15% error bars caused by fluctuations in the 
laser power and gas jet density. We extract the Gaussian beam parameters 
by fitting the measured intensity profile for each harmonic. 
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Figure 1: (a) Experimentally reconstructed spectral wavefront phase and amplitude 
at the slit position for different harmonics. The harmonics H15 (data points plotted 
as V) and H21 (data points plotted as V) are shown in red and blue respectively. 
The phase data points are interpolated by the dashed line, and the measured am- 
plitude data points are shown interpolated by solid line. The background signal 
arises from stray light, (b) The difference between the fitted parabolic phase and 
the measured phase. 



The wavefront phase <Pn(x, Zsiit)i where z s m is the slit position along 
the beam propagation and x is the coordinate across the slit, is parabolic 
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for all harmonic orders N (Fig. [T] (a)). Because the noise in the phase 
measurements is small, the error bars are not visible on the scale of Figure 
[H (a). To better estimate the phase error we subtract the fitted parabolic 
phase from the measured phase as shown in Figure [T] (b). Across the central 
part of the beam, containing 95% of the beam energy, the standard deviation 
of measured phase from ideal parabolic phase is 0.12rad. (better than A/50) 
for the 15th harmonic and 0.22rad (better than A/28) for the 21st harmonic. 
Details about the fitting procedure for amplitude and phase of the wavefront 
are provided in the Methods section. 



a 40 



N 30 

a 

.2 20 



o 
ci- 



10 







s 

o-10 
Pi 



(a) 


- - -*- - * - 


A" 


* 

* 

/ 5Z 






/ 

/ Ska 










1- 









13 15 17 1< 



Hamionic Order N 



21 23 25 



16 
14 
12 

bo 

8 
6 
4 



(b 



13 15 17 1< 



21 23 25 



hhg 



Harmonic Order N 



hhg 



Figure 2: (a) Relative waist position (SZ) as a function of harmonic order. The 
error bar is shown within 95% confidence bounds. The insert illustrates that dif- 
ferent harmonics "originate" from different positions in space. The meaning of the 
parameters SZ and Z\s is illustrated in the inset, (b) The beam waist size wqn as a 
function of harmonic order. The error bar is shown within 95% confidence bounds. 



Given the amplitude and the phase of the electric field of each harmonic 
at the scanning slit plane, we propagate the complex electromagnetic field 
backwards in time. Gaussian beams with parabolic wavefront phase remain 
Gaussian as they propagate. They can be characterized by their beam waist 
position Zn (relatively to the scanning slit) and their beam waist size ^on- 
Figure [2] (a) shows the reconstructed waist positions as a function of the 
harmonic order. The differences in the position Z-^ for different harmonics 
quantitatively show the different curvatures for each harmonic. The relative 
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waist position 5Z (Fig. [2] (a)) varies by more than their Rayleigh range 
{Zqn) which monotonically decrease from 11.8mm for H13 to 2.7mm for 

2 

H25 (Zotv = x° N ; where Atv is the wavelength of N th harmonic). The 
corresponding beam (Fig. [2](b)) waists vary from 15.2/^m for H13 to 5.3/um 
for H25. 
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Figure 3: Spatial-temporal coupling. Panel (e) shows the selected harmonics near 
their focal position assuming a focusing mirror of focal length f=100mm placed a 
the position of the scanning slit. We see that different harmonics focus at different 
positions and have different divergence. Panel (f) shows the Gouy phase of the 
same harmonics around the focus. Assuming transform limited pulse at the center 
of the focusing mirror = 0, R = - panel (d)) we can determine the pulses 
temporal profile at any position. Panels (a), (b) and (c) show respective pulses 
at Z — 155mm, R — (panel (a)); Z — 165mm, R — (panel (b)) and Z = 0, 
R = 300/i (panel (c)). In (a)-(d) the field strength is normalized to unity. Red, 
green and blue colors show the results for 13, 19 and 23 harmonic respectively. 



The dependence of wavefront structure on harmonic order means that it 
is impossible to focus the attosecond pulse into one diffraction limited spot 



6 



by achromatic optics (by using mirrors for example) . Figure [3] shows how 
the harmonics would focus if an f= 100mm mirror was placed at the scanning 
slit position (about 250mm from the jet). 

With different harmonics focusing at different longitudinal positions it 
is impossible to maintain a Fourier-limited attosecond pulses as the beam 
propagates. As soon as a particular harmonic component passes through 
the focus - it acquires a ir phase shift (the Gouy phase) . Figure [3] (f ) shows 
the Gouy phase across the focal region for different harmonics. 

Spectrally resolved wavefronts open the door for complete spatial-temporal 
reconstruction of attosecond pulses. Given the measured spectral wavefront 
(spectral wavefront phase - (j)^(r,z), and spectral wavefront amplitude - 
\E(f, z,ojn)\, where r - transverse coordinate across the beam, z - is the co- 
ordinate along the beam propagation) at any plane z = zq along the beam 
propagation, it is enough to have temporal information at any single point 
f = r~o across the beam to determine the temporal structure everywhere. 

Since it is sufficient to measure the spectral phase 9?jv(rb,zo) for any 
(fo,Zo), we could select the point, for example, by centering a small aper- 
ture on high-harmonic beam and measuring the spectral phase of the trans- 
mitted radiation using RABBIT, FROG-CRAB [9] or any other suitable 
technique. Then, the complete spatial (3D)+ temporal (ID) field profile of 
the beam can be reconstructed - 

E hhg (r,z,t) = ^2 \ E (r,z,u} N )\ exp i[<p N (r , z ) + <j>N{r,z) + Nu t] (1) 

where ujq is the angular frequency of the fundamental beam and the 
summation is over all emitted harmonic orders. Note, while the spectral 
wavefronts \E(r, z, ojn)\ 6xpi[(^jv(^^)] is measured across one plane (z = 
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zq), it can be found across any plane along the beam propagation using 
Fresnel-Kirchhoff diffraction formula [10] for each harmonic. 

To demonstrate this principle, we arbitrary assume that the attosecond 
pulses within the train are transform limited at the center of the focusing 
mirror (Fig. [3](d)) i.e. ip^ = Nuot e , where t e is harmonic emission time. All 
other parameters (i.e. the wavefront structure of each harmonic) are taken 
from our measurement. As we move along or across the beam, the relative 
spectral amplitude and phase of the harmonics changes. Consequently the 
temporal structure of each attosecond pulse in the pulse train depends on 
where the pulse is measured. Using the measured spectral wavefronts in 
our experiment, the temporal profiles is shown for selected points along the 
beam in Fig. E](a)-(c). Even in cases where the atto-chirp I12j has been 
corrected, it can only be corrected strictly speaking at one point. 

There are at least two approaches for bringing the harmonics to a com- 
mon focus and generating Fourier limited attosecond pulses and pulse trains. 
(l)The phase of any harmonic depends on both the intensity and the phase 
of generating beam. Optimizing the spatial profile (both phases and am- 
plitudes) of the fundamental beam will allow us to minimize the chromatic 
wavefront variations between different harmonics. (2)The chromatic aberra- 
tions that we measure are reproducible. They can, in principle, be corrected 
with specially designed linear optics. 

In conclusion, within the signal-to-noise limits of our experiment, we have 
obtained complete (phase and amplitude) information on wavefront profile 
of each harmonic. We have shown that spectrally resolved wavefronts allow 
the field distribution to be found anywhere. If we know the position of 
the generating medium with sufficient accuracy, this includes the position 
at the generating medium itself. The wave front amplitude and phase in 
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the generating medium contains an imprint of the underlying single atom 
or molecule response. Within a 3-step model [13] , for a thin, low density 
gas, the measured wavefront phase (j)]sf{f,Zj e t) is given by the sum of four 
contributions: 

<j) N (r, Zjet) = 4>t{I{r, z jet ),u) N )-\a N \I(r, z jet )+(J) r (u N , z je t)+N(j)f und (f, z jet ) 

(2) 

Where 4>{ UQ d(r, Zj e t) is the wavefront phase of the fundamental beam at the 
jet position. (f)f un <i(r, Zj et ) can be measured. 4> r (uj]y, Zjet) is the phase contri- 
bution of the transition moment [14]. <ft r (ujN, zj e t) is intensity independent 
for a single orbital, but intensity dependent for multiple orbitals. The phase 
that the electron acquires in the continuum is —\aj\r\I(f,Zj e t). This phase 
can be calculated within the strong field approximation and its corrections 
[3 E]. (j)t(I(f, Zjet),uJN) is the yet not fully understood tunneling phase 
I16j . It might be intensity dependent. Therefore, wavefront structure 
will serve as a sensitive measurement of the intensity dependence of these 
parameters. 

Once the relative phase of each harmonic is determined, then we have 
shown that we know the temporal profile of an attosecond pulse anywhere 
in space. No matter what the pulse temporal structure should prove to 
be, we have already determined that there will be significantly different 
temporal profile at the center and at the edges of a beam and the temporal 
profile must change as the beam passes through the respective foci of the 
harmonics. Therefore, any spatially extended measurement with attosecond 
pulses will be affected and any measurement of the pulse itself is also affected. 
The result of any experiment will be effectively averaged over the different 
temporal profiles at different positions in the interaction region of the focal 
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volume. 

Finally, attosecond science has been restricted to spatially averaged mea- 
surements to image orbitals [T7], tracing molecular dynamics [18], identify 
and time resolve tunneling wave packets [15] and to follow Auger decay [19] . 
Spectrally resolved wavefronts and the complete spatio-temporal character- 
ization of attosecond pulses that they facilitate will allow us much greater 
experimental precision in all of these experiments. 

We thank Mike Spanner, Misha Ivanov, Kyung Teac Kim and Eric 
Constant for stimulating discussions. We acknowledge financial support 
of Canada's National Research Council, National Sciences, Engineering Re- 
search Council, the Canada Research Chair Program and MURI grant W911NF- 
07-1-0475. E.F acknowledges the support of Marie Curie International Out- 
going Fellowship and H.N acknowledges the support of JST's Presto pro- 
gram. 

0.1 Methods 

The laser beam from Ti:Sapphire amplifier was spatially filtered by propa- 
gating it through hollow-core fiber in vacuum and the beam was then focused 
with an //# = 80 lens onto the pulsed gas jet. 

The gas jet in which the harmonics were generated has 250/im aperture 
and was operated with a backing pressure of 2.7atm. The harmonic signal 
was recorded using a 40mm diameter Burle imaging microchannel plate 
(MCP) in the imaging plane of the spectrometer and the back-side phosphor 
screen was imaged onto the CCD camera. Further details on experimental 
setup can be found in Ref. [8j. 

The Rayleigh range of the fundamental beam (Zr ~ 10mm) was much 
longer than the interaction region of the jet (d ~ 0.4mm) where harmonics 
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were generated. The intensity of the fundamental beam was kept moderate 
to prevent excessive ionization and plasma creation. This configuration is 
particularly useful when we aim to study single atom/molecular response 
and has been widely used in attosecond science [17] . 

The high harmonic radiation was diffracted through a horizontal scan- 
ning slit positioned in front of the XUV spectrometer [8]. The spectrometer 
is built with variable groove spacing flat field Hitachi imaging grating (model 
001-0266). The system resolves the spectrum in horizontal direction, and 
allows essentially free-space field propagation in vertical direction. Two- 
dimensional images are taken for each position of the scanning slit using 
imaging multichannel plate (MCP) and cooled CCD camera. 

The relative vertical position of the diffraction pattern's centroid is pro- 
portional to the wavefront slope of the sampled wavefront slice. We re- 
construct the wavefront phase in Fig. [1] (a) using the following relation: 
5j = Sz ■ d Zj , where yj - is the relative position of the centroid of the 
diffracted pattern, Zj - relative slit position, d - is the distance between the 
scanning slit and the imaging plane, Sz - is the scanning slit iteration step, 
and Sj - is the optical path difference across the slit. Then we determine 
the sampled wavefront phase profile - 4>N exp for the harmonic wavelength - 
Xn across the scanning direction using <fi l Nexp = j^- J2)=o We determine 
the wavefront amplitude P exp at each sampling point by integrating inten- 
sity of the corresponding diffraction pattern. Further details of SWORD are 
discussed in Ref. [HJ. 

The experimental wavefront phase {(j)% exp } sampled at positions {x l exp } 
was fit with parabolic function 4>n{x) = aN(x — xo) 2 for each harmonic. The 
Nonlinear Least Square method was used for the fit. To minimize the noise 
contribution, the phase is fit within the range where measured intensity 
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is larger than 5% relative to the peak intensity. The wavefront intensity 
distribution {P exp \ measured at the scanning slit was fit to the Gaussian 
profile I(x) = Iq exp[— 2{x — xq) 2 /w 2 sN ] + J& , where 1$ is maximum intensity, 
w s n - Gaussian width at the slit position, and lb - background signal. The 
measured parabola curvature - a^v, and the Gaussian width - w s n for each 
harmonic (H13-H25) were used to retrieve the virtual waist position Z^ and 
waist size wqn. A sensitivity analysis showed that the Zn is highly sensitive 
to cijv and depends weekly on w s n (ex.: 10% change in the Gaussian width 
causes only ~ 0.1% change in Zjy). We use the confidence bound of 95% 
in determination of the error bars for the relative waist position 8Z with 
respect to the scanning slit as shown in Fig [2 

We have measured short trajectory harmonics. The divergence of the 
short trajectories mainly depends on the position of the jet relative to focal 
position of the fundamental beam. For all jet positions we find that the 
virtual origin of the harmonics differs by more than their confocal parameter 
across the spectrum. 
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